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Abstract 

Tuning the finite element model using measured 
data to minimize the model uncertainties is a 
challenging task in the area of structural 
dynamics. A test validated finite element model 
can provide a reliable flutter analysis to define 
the flutter placard speed to which the aircraft 
can be flown prior to flight flutter testing. 
Minimizing the difference between numerical 
and experimental results is a type of 
optimization problem. Through the use of the 
National Aeronautics and Space Administration 
Dryden Flight Research Center’s (Edwards, 
California, USA) multidisciplinary design, 
analysis, and optimization tool to optimize the 
objective function and constraints; the mass 
properties, the natural frequencies, and the 
mode shapes are matched to the target data and 
the mass matrix orthogonality is retained. The 
approach in this study has been applied to 
minimize the model uncertainties for the 
structural dynamic model of the aerostructures 
test wing, which was designed, built, and tested 
at the National Aeronautics and Space 
Administration Dryden Flight Research Center. 
A 25-percent change in flutter speed has been 
shown after reducing the uncertainties . 

1 Introduction 

A test article called the aerostructures test wing 
[ATW] was developed and flown at the 
National Aeronautics and Space Administration 
[NASA] Dryden Flight Research Center 
[DFRC] (Edwards, California, USA) on the 
NF15B (McDonnell Douglas, now The Boeing 
Company, Chicago, Illinois, USA) test bed 
aircraft, as shown in figure 1 , for the purpose of 


demonstrating and validating flutter prediction 
methods during flight [1], The first 
aerostructures test wing [ATW1], flown in 
2001, was originally developed to directly 
address requests for better flight flutter test 
techniques by providing a functional flight test 
platform. While the first series of tests was 
extremely successful, the minimum amount of 
instrumentation (structural accelerometers and 
strain gages) was chosen to satisfy the scope of 
the program. These sensors were limited in the 
capability to answer questions of aeroelastic 
interactions, sources of nonlinearity, physical 
mechanisms of aeroelastic coupling, and 
feedback dynamics between the structure and 
aerodynamics. 

A second aerostructures test wing 
[ATW2], as shown in figure 2, has been built 
for the demonstration of state-of-the-art sensor 
technologies for simultaneous, distributed, and 
collocated measurement of shear stress (skin 
friction); steady and unsteady pressures; and 
structural strain and accelerations for mode 
shapes and other modal properties. This wing 
was flown on the NF15B in December 2009. 

A block diagram for a robust flutter 
analysis procedure used at NASA DFRC is 
given in figure 3. Using a finite element [FE] 
model for a structural dynamic analysis 
becomes increasingly important in the modem 
aircraft design and analysis processes. In 
general, the quality of an initial FE model of an 
aircraft is not guaranteed, so we need to perform 
the ground vibration test [GVT] to validate the 
FE model. In most cases, these newly built FE 
models need to be tuned to minimize 
uncertainties in the structural dynamic FE 
models and flutter boundary results. Robust 
flutter analyses performed at NASA DFRC are 
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mainly based on these validated FE models as 
shown in figure 3. 

Model tuning is a common method to 
improve the correlation between numerical and 
experimental modal data, and many techniques 
have been proposed [2], Literature reviews on 
finite element model tuning are summarized in 
reference [3]. These techniques can be divided 
into two categories: direct methods (adjust the 
mass and stiffness matrices directly) and 
parametric methods (correct the models by 
changing the structural parameters). The direct 
methods correct mass and stiffness matrices 
without taking into account the physical 
characteristics of the structures and may not be 
appropriate for use in model tuning processes. 
In this paper, the tuning method used in the 
optimization process is the parametric method. 
In the optimization process, structural 
parameters are selected as design variables: 
structural sizing information (thickness, cross 
sectional area, area moment of inertia, torsional 
constant, et cetera), point properties (lumped 
mass, spring constants, et cetera) and materials 
properties (density, Young's modulus, et cetera). 
Objective function and constraint equations 
include mass properties, mass matrix 
orthogonality, frequencies, and mode shapes. 
The use of these equations minimizes the 
difference between numerical results and target 
data. 

The primary objective of this study is to 
reduce uncertainties in the structural dynamic 
FE model of the ATW2 to increase the accuracy 
of flutter speed prediction. Discrepancies are 
common between the test data and numerical 
results. However, the FE model can be fine 
tuned through the use of GVT data. Accurate 
and reliable GVT results are important to this 
adjusting process. 

The secondary objective of the current 
study is to add model-tuning capabilities [3] in 
NASA Dryden’s object-oriented 

multidisciplinary design, analysis and 
optimization [MDAO] tool [4] . This model 
tuning technique [3] is essentially based on a 
non-linear optimization problem. 


2 Object-Oriented MDAO Framework 

The heart of the object-oriented MDAO tool is 
the central executive module [CEM] as shown 
in figure 4. In this module the user will choose 
an optimization methodology, provide side 
constraints for continuous as well as discrete 
design variables and external file names for 
performance indices; which communicate 
between the CEM and each analysis module, 
submit script commands to prepare input data 
for each analysis code, execute analyses codes, 
compute performance indices using post- 
processor codes, and compute an objective 
function and constraints values from 
performance indices. The CEM was written in 
FORTRAN, and script commands for each 
performance index were submitted through the 
use of the FORTRAN call system command. 

The performance indices for structural 
optimization problems can include total weight, 
safety factors, natural frequencies, mode shapes, 
flutter speed, divergence speed, structural 
responses at sensor locations, et cetera. On the 
other hand, the performance indices for the FE 
model tuning problems can include off-diagonal 
terms of orthonormalized mass and stiffness 
matrices, the total error between the measured 
and computed mode shapes at given sensor 
points, errors between computed and measured 
frequencies, total weight, center of gravity [CG] 
locations, moment of inertia, et cetera. 

Two optimization codes are included in 
the object-oriented MDAO tool: design 

optimization tools [DOTs] [5] based on a 
gradient-based algorithm [6] and a genetic 
algorithm [GA] [7], The MSC/NASTRAN 
(MSC Software Corporation, Santa Ana, 
California, USA), for example, uses a gradient- 
based approach for optimization [8] . A 
drawback to this approach is the necessity to 
compute finite difference or analytical 
sensitivity values to perform the search, which 
often requires prior experience based on input 
defining the problem and search directions. 

The DOT is a commercial optimization 
code that can be used to solve a wide variety of 
nonlinear optimization problems. When the 
DOT requires the values of the objective and 
constraint functions corresponding to a 
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proposed design, it returns control to the CEM. 
The CEM calls the DOT again to obtain the next 
design point. This process is repeated until the 
DOT returns a parameter to indicate that the 
optimum objective function is reached. 
Gradient-based algorithms work well for 
continuous design variable problems, whereas, 
GAs can easily handle continuous as well as 
discrete design variable problems. When there 
are multiple local minima, GAs are able to find 
the global optimum results, whereas, gradient- 
based methods may converge to a locally 
minimum value. 

The GA is directly applicable only to 
unconstrained optimization so it is necessary to 
use some additional methods to solve the 
constrained optimization problem. The most 
popular approach is to add penalty functions, in 
proportion to the magnitude of the constraint 
violation, to the objective function [9]. 

The general form of the exterior penalty 
function is shown in equation (1): 

<7 _ r 

L(x) = F(x)+ £ A,g, (x) + £ V< 7 h j ( X ) (1) 

(=1 j=\ 

where l(x) indicates the new objective 
function to be optimized, F(x) is the original 
objective function, g,- (x) is the inequality 
constraint, h ; (x) is the equality constraint, A,- 

are the scaling factors, X is the design variables 
vector, and q and r are the number of inequality 
and equality constraints, respectively. 

Five optimization methodologies 
available in the CEM are described as follows: 

• Genetic algorithm with continuous, 
discrete or mixed (continuous/discrete) 
design variables [GCD] ; 

• Gradient-based algorithm, that is DOT, 
with continuous design variables [DC] ; 

• Start with GCD then continue with DC; 

• Start with GCD then continue with DC 
then continue with genetic algorithm 
with discrete design variables [GD]; and 

• Start with DC then continue with GD. 


3 Structural Dynamic Model Tuning 
Procedure 

Discrepancies in frequencies and mode shapes 
are minimized using a series of optimization 
procedures [3, 10, 11]. The numerical mass 
properties, the mass matrix orthogonality, and 
the natural frequencies and mode shapes are 
matched to the target values based on the 
following three tuning steps. 

3.1 Step 1: Tuning Mass Properties 

The difference in the numerical and target 
values of the total mass, the center of gravity 
[CG] location, and mass moment of inertias at 
the CG location are used as performance 
indices. The resulting 10 performance indices, 
as shown in equations (2) through (11), are 
defined to minimize the uncertainties in the 
rigid body dynamics. 


J 1= (W-W G ) 2 /W 2 (Total weight) 

(2) 

j 2 =(x-x g ) 2 /x 2 

(X-CG) 

(3) 

J 3 =(y-y g ) 2 /y 2 

(Y-CG) 

(4) 

j 4 =(z-z g ) 2 /z 2 

(Z-CG) 

(5) 

J 5 = ( X XX “ X XXG ) 2 1 X XXG 

(ixx at CG location) 

(6) 

J 6 = ( X YY _I YYg) 2 /X YYG 

(iyy at CG location) 

(7) 

J 7 =( I ZZ _I ZZ G ) 2 /X ZZG 

(izz at CG location) 

(8) 

J 8 =( I XY _I XYg) 2 /X XYG 

(ixy at CG location) 

(9) 

J 9 =( I YZ _I YZg) 2/I yZG 

(iyz at CG location) 

GO) 

j io = ( x zx _i zxg) 2 /X ZXG 

(izx at CG location) 

GD 


In the case of a DOT optimization, it is 
recommended to start an optimization procedure 
at a feasible domain. When the optimization 
procedure starts in an infeasible domain, there is 
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no guarantee that the design will move to a 
feasible domain. To this end, these ten sub- 
optimization criteria are used (if test data is 
available) to make a feasible starting 
configuration for the second and the third 
optimization steps. 

The exterior penalty function approach 
given in equation (1) is used for the GA 
optimization. For example, the performance 
index and performance indices J 2 through 
J |0 can be selected as an objective function and 
constraints, respectively. Therefore, the 
following objective function, J , can be used for 
the GA optimization as shown in equation (12): 

10 

J = J\ + ^ ^i'li (12) 

(=2 

where A,- is the scaling factors and i = 2,. . . , 10. 


3.2 Step 2: Tuning Mass Matrix 

The off-diagonal terms of the orthonormalized 
mass matrix are reduced to improve the mass 
orthogonality as shown in equation (13): 


and difficulties in measuring the rotational 
DOFs. 


3.3 Step 3: Tuning Frequencies and Mode 
Shapes 

Two options can be used for tuning the 
frequencies and mode shapes. In the first option, 
the objective function considered combines a 
performance index J| 2 (the normalized errors 
between GVT and computed frequencies) as 
shown in equation (15) with a performance 
index J 13 (the total error associated with the 
off-diagonal terms of the orthonormalized 
stiffness matrix) as shown in equation (16). 


J 12 = E 

i = 1 


Q/ -Q/G f 

£2 ( - J 


(15) 


J 13= E ( K ij ) (option 1) (16) 

i=l,j=l,i*j 

The matrix K is obtained from the 
following matrix products as shown in equation 
(17): 


Jll= E ( M ‘j) (13) 

i=Uj=U*j 


K = 0^r r KTd) G 


(17) 


where n is the number of modes to be matched 
and M is defined as shown in equation (14). 

M = O^T 7 ’MT(l) G (14) 

The performance index J n is used as 
the objective function, and performance indices 
J| through J 10 are used as constraints in this 
step. In equation (14), the mass matrix M is 
calculated from the FE model while the target 
eigen-matrix <t> G is measured from the GVT. 
The eigen-matrix <t> G remains constant during 
the optimization procedure. A transformation 
matrix T [3] in equation (14) is based on Guyan 
reduction, improved reduction system [12] or 
system equivalent reduction expansion process 
[SEREP] [13]. This reduction is required due to 
the limited number of available sensor locations 


where the stiffness matrix, K , is calculated 
from the FE model. 

In the second option, the objective 
function considered combines the performance 
index J| 2 with a new performance index J J3 
(the total error between GVT and computed 
mode shapes at given sensor points) as shown in 
equation 18. 

m n 2 

J 13 = E EK-^'/g) (option 2) (18) 

i= 1 7=1 

The performance index J 13 , in the 
second option, is much simpler than in the first 
option for this application. Any errors in both 
the modal frequencies and the mode shapes are 
minimized by including an index for each of 
these in the objective function. For the second 
option, a small number of sensor locations can 
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be used at which errors between the GVT and 
computed mode shapes are obtained. 

Instead of using the summation in 
equations (13) and (15), each individual term 
can be a separate performance index. In this 
case, the total number of performance indices 
for steps 2 and 3 will increase, however, it is 
easier to apply the military standard 

requirements. 

Any one or combinations of 

performance indices J| through J 13 can be used 
as the objective function with the other 
performance indices treated as constraints. This 
gives the flexibility to achieve the particular 
optimization goal while maintaining the other 
properties as close to the desired target value as 
possible. The optimization problem statement 
can be written as 

Minimize J, 

Such that I* - e k , for k= 1 through 13 and k * i 

where e k is a small value that can be adjusted 
according to the tolerance of each constraint 
condition. 

4 Test Article 

The ATW2 was used to demonstrate NASA 
Dryden’s object-oriented MDAO tool through 
the process of ground vibration testing and the 
model tuning technique. This test article was a 
small-scale airplane wing comprised of an 
airfoil and wing tip boom as shown in figure 5. 
The major structural component was the 
composite wing spar with a rectangular cross 
section. This wing had equally spaced four ribs 
and composite upper and lower cover skins. The 
empty space between the upper and the lower 
cover skins was filled with plastic foam. This 
wing was based on a NACA-65A004 airfoil 
shape with a 3.28 aspect ratio. The wing had a 
half span of 18 in (45.720 cm) with a root chord 
length of 13.2 in (33.528 cm) and a tip chord 
length of 8.7 in (22.098 cm). The total area of 
this wing was 197 in 2 (1270.965 cm 2 ). The wing 
tip boom was a 1-in (2.540 cm) diameter hollow 
tube of 21.5 in (54.610 cm) in length. The total 
weight of the wing was 2.66 lb (0.9928 kg). 


Since the ATW2 was attached to the 
NF15B flight test fixture, the construction of the 
wing was limited to lightweight materials with 
no metal for the safety of re-contact with the 
aircraft after a possible separation. The wing 
and spar were constructed from fiberglass cloth, 
the tip boom used carbon fiber composite, the 
wing core consisted of rigid foam, and the 
components were attached by epoxy. The wing 
skin was made of three plies of fiberglass cloth 
with a 0.01 -in (0.254 mm) thickness. The 
internal spar located at the 30-percent chord line 
was composed of 10 plies with a 0.05-in 
(1.270-mm) thickness of carbon at the root, and 
decreasing to 1 ply 0.005-in (0.127 mm) 
thickness at the tip. 

5 Test Setup 

Ground vibration tests were performed to 
determine the dynamic modal characteristics of 
the ATW2. In the test set up, the ATW2 was 
clamped on to a circular aluminum plate, which 
was bolted to a mounting panel, and then 
installed into a small strong back called the 
ground test fixture in the NASA Dryden Flight 
Loads Laboratory. The PONTOS 
photogrammetry optical measuring system 
(Gesellschaft fur Optische Messtechnik, 
Braunschweig, Germany) [14], as shown in 
figure 6, was used to measure output 
displacement at the sensor points, shown in 
figure 7. For the excitation method, an impulse 
hammer with an impedance head was used to 
excite the ATW2’s natural frequencies and 
mode shapes as well as to measure input forces. 

PONTOS is a non-contact optical 3D 
measuring system. It analyzes, computes, and 
documents object deformations, rigid body 
movements, and the dynamic behavior of a 
measuring point [14]. The PONTOS system 
provides an alternative for complex sensor 
technology like laser sensor, draw- wire sensors, 
or accelerometers, which are commonly used in 
the GVT for measuring responses of the 
structure. The features of the PONTOS system 
include: 

• an unlimited number of sensors (The 
sensor markers are extremely 
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lightweight so that a large number of 
sensors can be used at the same time 
without essentially altering the total 
weight or the mode shapes of the 
structure.), 

• non-contact acquisition of the precise 3D 
position of any number of measuring 
points, 

• mobility and flexibility due to an easy 
and compact measuring system, and 

• easy and quick adaptation to different 
measuring volumes and measuring tasks. 

The limitations of the PONTOS system include: 

• measuring frame rate up to 500 Hz at 
1280x1024 pixels, 

• measuring volume up to 1700 by 1360 
by 1360 mm 3 , and use on a plane or a 
slightly curved surface. 

6 Flutter Analysis Before Model Tuning 

Natural frequencies obtained from the GVT and 
the MSC/NASTRAN modal analyses using FE 
models before model tuning are compared in 
table 1. The starting configuration of the FE 
model was obtained from the ATW1 model. 
The GVT frequencies in this table were based 
on the time history responses data, as shown in 
figure 8, collected by the PONTOS system at 
each of the sensor points. The eigensystem 
realization algorithm routine, which was 
developed by Juang and Pappa [15] at NASA 
Langley Research Center, was then used to 
identify the frequencies and mode shapes of the 
system. The FE model for the MSC/NASTRAN 
[16] modal analysis is shown in figure 9. 
Detailed GVT sensor locations and GVT mode 
shapes are provided in reference [17]. 
Corresponding natural frequencies and mode 
shapes computed using the MSC/NASTRAN 
code are shown in figure 10. In table 1 and 
figure 10, the mode numbers 1,2, and 3 are the 
first bending, first torsion, and second bending 
modes, respectively. 

A matched flutter analysis of the ATW2 
using the ZAERO code [18] is presented in this 
section. The first 10 natural modes were used 
for the flutter analysis. The aerodynamic model 
has 200 surface elements on the wing and 216 


body elements on the boom. Forty-four splining 
points were used between the structural 
dynamic and aerodynamic models. The 
aerodynamic model of the ATW2 and the first 
three splined mode shapes on the aerodynamic 
model are shown in figures 11 and 12, 
respectively. The matched flutter analyses were 
performed at Mach 0.60, 0.75, 0.82 and 0.95. 
The aerodynamic influence coefficient matrices 
at each Mach number were generated at 16 
reduced frequencies and the g-method [18] was 
used in the matched flutter analysis. 

The speed versus damping, V-g, and 
speed versus frequency, V-w, curves from the 
matched flutter analysis at Mach 0.82 (the 
ATW1 flutter Mach number) [1] before model 
tuning are given in figure 13. The structural 
damping used for the flutter speed computation 
was 3 percent. The flutter speed was 

407.4 KEAS, flutter frequency was 22.86 Hz, 
and the corresponding altitude was 15,010 ft 
(4,575 m), as shown in table 2. 

Modal participation factors for the first 
flutter mode are given in table 3 . Modal 
participation of the first three natural 

frequencies was more than 96 percent before 
model tuning. In-plane modes, mode numbers 4, 
6, and 8, did not contribute to the first flutter 
mode. In table 3, we can conclude that the 
primary structural dynamic modes for the first 
flutter mode were modes 1, 2 and 3 and the 
secondary modes for the first flutter mode were 
modes 5, 7, 9 and 10. The first mode is the 
dominant mode for the primary flutter mode. 

7 Model Tuning 

Based on military standards [19, 20], the 
frequency error should be less than 3 percent for 
the primary modes and 10 percent for the 
secondary modes. Using a norm of frequency 
differences as an objective function with mass 
properties, mass orthogonality, and mode shapes 
as constraint equations, the frequencies after 
model tuning are presented in table 1. Before 
model tuning the frequency error for the second 
mode was 47.3 percent and the frequency error 
for the third mode was 11.9 percent. These 
frequency errors violate the 3-percent frequency 
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error requirements in military standards. After 
model tuning the 2.09-percent error becomes 
3.19 percent, the 47.3-percent error becomes 
1.38 percent, and the 11.9-percent error 
becomes 0.39 percent. After model tuning, 
3 percent frequency error requirements for the 
primary modes were satisfied for modes 2 and 
3 , and close to being satisfied for mode 1 . 

Since Guyan reduction is a static 
condensation, it is only accurate for lower 
modes. For higher modes, the errors become too 
large as shown in table 1 . Unlike Guyan 
reduction, the SEREP preserves the dynamic 
character of the original full system model for 
selected modes of interest. The dynamic 
characteristics of the reduced model were 
virtually the same as the full model as shown in 
table 1. Therefore, the SEREP model reduction 
process was chosen for this ATW2 model 
update application. 

Table 4 shows the total weight, 
orthonormalized mass matrix, and modal 
assurance criterion [MAC] [3] values of the 
ATW2 before and after model tuning. Based on 
military standards, the off-diagonal terms of the 
orthonormalized mass matrix should be less 
than 10 percent [19, 20]. The off-diagonal terms 
of the orthonormalized mass matrix, with a 
maximum of 38 percent before model tuning, 
are minimized in the second tuning step. The 
maximum off-diagonal term of 7.43 percent 
after model tuning is observed in table 4 and 
satisfies the 10-percent limitation in military 
standards. Model correlation with test data prior 
to model tuning was poor and unacceptable to 
proceed to flight. The MAC values of 0.70 and 
0.75 for modes 2 and 3 before model tuning 
become 0.97 and 0.95, respectively. Therefore, 
we can conclude that excellent model 
correlation with the test data was achieved after 
model tuning, which leads to a more reliable 
flutter speed prediction. 

8 Flutter Analysis After Model Tuning 

Typical V-g and V-w curves from the matched 
flutter analysis at Mach 0.82 after model tuning 
are given in figure 14. The same 3-percent 
structural-damping value is used for flutter 


speed computation. In figure 14 we can observe 
a steeper merging behavior for the first bending 
and torsion modes after model tuning. It is 
important for the safety of flight that the model 
tuning results predict a more sudden, rather than 
a more gradual, flutter onset. 

The results of the matched flutter 
analysis before and after model tuning are 
summarized in tables 2 and 5. The modal 
participation factors in table 5 indicate that the 
first three modes have more than a 95 -percent 
contribution to the first flutter mode before 
model tuning and more than a 98-percent 
contribution after model tuning. The second 
mode, which has the maximum frequency 
difference before model tuning as shown in 
table 1, becomes more than 91 percent of the 
modal participation after model tuning at 
Mach 0.60. Therefore, the second mode after 
tuning becomes the dominant mode for the 
primary flutter mode. The resulting flutter speed 
difference due to model uncertainty was also a 
maximum at Mach 0.60, as shown in table 2. 
Table 2 shows that the flutter speed difference 
varies from 25.4 percent at Mach 0.60 to 
3.0 percent at Mach 0.95. As shown in tables 2 
and 5 after model tuning, the modal 
participation of the second mode, which is the 
dominant mode, decreased when Mach number 
increased; the differences in flutter speed due to 
uncertainty also decreased when Mach number 
increased. 

Finally, flutter boundaries of the ATW2, 
before and after model tuning, are compared 
with the flight envelopes as shown in figure 15. 
The solid and long dash dotted lines are the 
flight envelope of the mother ship and its 15- 
percent margin, respectively. The dashed line 
represents the ATW 2 test envelope that is 
planned for flight, and the long dash double dot 
line is the 15 -percent margin of the ATW2 test 
envelope. The solid line with the square marker 
is the flutter boundary before model tuning and 
the solid line with the diamond marker 
represents flutter boundary after model tuning. 
It should be noted that flutter boundary after 
model tuning is more conservative than before 
model tuning. 
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9 Concluding Remarks 

This paper describes the reduced uncertainty 
procedures for the ATW2, which was developed 
at the NASA Dryden Flight Research Center for 
demonstrating flutter and advanced aeroelastic 
test techniques. It has been shown that model 
uncertainties can be reduced through the tuning 
of the finite element model using the NASA 
DFRC’s object-oriented MDAO tool. 

After tuning the FE model, the 
frequency differences between GVT and 
analytical results are mostly within 3 percent 
and the off-diagonal terms of the 
orthonormalized mass matrix are within 

10 percent satisfying the military standards. 
Excellent mode shape correlations were also 
achieved through the high MAC value (greater 
than 95 percent). With the updated FE model, 
the flutter speed prediction can be reduced as 
much as 25 percent at Mach 0.60. 
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REDUCED UNCERTAINTIES IN THE FLUTTER ANALYSIS OF THE AEROSTRUCTURES TEST WING 


Tables 

Table 1. Natural frequencies of the ATW2 before and 
after model tuning. 


Before tuning After tuning 


Mode 

(Hz) 

MSC/NASTRAN 

Error 

MSC/NASTRAN 

Error 


(Guyan/full; Hz) 

(%) 

(SEREP/full; Hz) 

(%) 

1 

17.24 

17.61/17.60 

2.15/ 

2.09 

17.79/17.79 

3.19/ 

3.19 

2 

44.10 

23.27/23.26 

-47.20/ 

-47.30 

44.71/44.71 

1.38/ 

1.38 

3 

84.00 

99.02/93.99 

17.90/ 

11.90 

84.33/84.33 

0.39/ 

0.39 


Table 2. Flutter boundaries before and after model tuning. 


Mach Mach Mach Mach 

= 0.60 =0.75 =0.82 =0.95 



Speed, 

KEAS 

453.0 

421.5 

407.4 

377.9 

Before 

tuning 

Frequency, 

Hz 

23.18 

22.97 

22.86 

22.53 


Altitude, 

ft 

-7,501 

8,751 

15,010 

25,590 


Speed, 

KEAS 

361.3 

364.6 

365.3 

367.0 

After 

tuning 

Frequency, 

Hz 

38.47 

37.73 

37.32 

36.20 


Altitude, 

ft 

5,101 

16,080 

20,310 

26,920 

Speed difference, 
percent 

25.4 

15.6 

11.5 

3.0 


Table 3. Modal participation factors at Mach 0.82 before 
model tuning. 


Table 4. Summary of total weight, orthonormalized mass 
matrix, and MAC values for the ATW2 before and after 
model tuning. 



Measured 

Before tuning 

After tuning 

Total weight 

2.66 lb 

1 .76 lb (error 34%) 

2.72 lb (error 2.3%) 



1 2 3 

1 2 3 

Orthonormalized 

1.0000-0.2490 0.3800 

1.0000 0.0395-0.0565 

mass matrix 2 

-0.2490 1.0000-0.6610 

0.0395 1.0000-0.0743 


3 

0.3800-0.6610 1.0000 

-0.0565-0.0743 1.0000 

Modal 

Mode 1 

0.97 

0.99 

assurance 

Mode 2 

0.70 

0.97 

criteria 

Mode 3 

0.75 

0.95 


Table 5. Modal participation factors before and after 
model tuning. 


Mode 

Frequency, 

Modal participation factors before model tuning, 
percent 

Hz 

Mach Mach 

Mach 

Mach 



= 0.60 =0.75 

= 0.82 

= 0.95 

1 

17.60 

68.1 72.9 

75.0 

79.7 

2 

23.26 

22.2 95.5* 18.3 96.2* 16.8 96.6* 

13.6 97.6* 

3 

93.99 

5.2 5.0 

4.8 

4.3 

Mode 

Frequency, 

Modal participation factors after model tuning, 
percent 

Hz 

Mach Mach 

Mach 

Mach 



= 0.60 =0.75 

= 0.82 

= 0.95 

1 

17.79 

6.7 9.8 

12.3 

24.3 

2 

44.71 

91.9 99.2* 88.4 99.1* 85.699.0* 

72.4 98.3* 

3 

84.33 

0.6 0.9 

1.1 

1.6 


*: These are the sum of modes 1 , 2 and 3 for each Mach number. 


Mode 

Frequency, Hz 

Modal participation factor, 
percent 

1 

17.60 

75.0 


2 

23.26 

16.8 

96.6* 

3 

93.99 

4.8 


4 

135.40 

0.0 


5 

163.10 

2.6 


6 

174.50 

0.0 


7 

257.50 

0.5 

3.4** 

8 

391.60 

0.0 


9 

394.30 

0.1 


10 

445.60 

0.3 


*: The total 

of the first three 

modes 


**: The total 

of modes 5, 7, 9 

and 10 



Figures 



Fig. 1. The aero structures test wing mounted on the 
NF-15B for flight flutter testing. 
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Fig. 2. The aero structures test wing 2. 


Measure weight CG, 
moment of inertia, 
and GVT data 


Fig. 3. Robust flutter analysis procedure at NASA DFRC. 



Fig. 4. Central executive module. 



Finite element model 

■ 

i 

1 

MDAO tool with model 
tuning capability 


1 

Validated 

finite element model 


i 

Create unsteady 
aerodynamic model 


i 

Perform flutter analysis 

1 

1 



Fig. 6. The PONTOS photo grammetry optical measuring 
system. 



Fig. 7. GVT sensor and excitation locations. 
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REDUCED UNCERTAINTIES IN THE FLUTTER ANALYSIS OF THE AEROSTRUCTURES TEST WING 


a> 

CO 

s 

o 

CL 

CO 

0 ) 


Rec. 13: 116X+ 112Y+ 



0 ) 

8 180 
ro . 

£ -180 


Rec. 13: 116X+ 112Y+ 



■§ .001 
Uzi 


100 150 

Frequency, Hz 


250 


Fig. 8. Typical time history and frequency response GVT 
results for the ATW2. 




Fig. 10a. Mode 1 (1 st bending): 17.60 Hz. 



Fig. 10b. Mode 2 (1 st torsion): 23.26 Hz. 



Fig. 10c. Mode 3 (2 nd bending): 93.99 Hz. 
Fig. 10. Full model mode shapes before tuning. 



Fig. 11. ZAERO unsteady aerodynamic model. 



Fig. 12a. Splined mode 1. 




Fig. 12c. Splined mode 3. 


Fig. 12. Splined mode shapes on unsteady aerodynamic 
model. 


11 



CHAN-GI PAK AND SHUN-FAT LUNG 



200 300 400 500 600 700 

Speed, KEAS 

Fig. 13a. V-g plots. 
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Fig. 14b. V-w plots. 
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Fig. 13b. V-w plots. 
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700 


Fig. 13. V-g and V-w plots for the ATW2 at Mach 0.82 
before model tuning (plot the first 5 modes) . 



Speed, KEAS 

Fig. 14a. V-g plots. 


Fig. 14. V-g and V-w plots for the ATW2 at Mach 0.82 
after model tuning (plot the first 5 modes) . 

Mother ship envelope 

Mother ship envelope x 1.15 

ATW2 test envelope 

ATW2 test envelope x 1.15 

■ Flutter boundary before model tuning 
Flutter boundary after model tuning 



Fig. 15. Flutter boundary before and after model tuning. 
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Reduced Uncertainties in the Flutter Analysis 
of the Aerostructures Test Wing 

Chan-gi Pak, Ph.D* and Shun-fat Lung, Ph.D** 

*: NASA Dryden Flight Research Center 
**: Contractor, Tybrin Corporation 



Objectives 


□ Support the Aeronautics Research Mission Directorate (ARMD) guidelines at NASA's Dryden 
Flight Research Center. 

❖ Supported by Subsonic Fixed Wing (SFW) & Supersonics (SUP) projects under Fundamental Aeronautics 
(FA) Program 

□ Reduce uncertainties in the structural dynamics model of the Aerostructures Test Wing 2 to increase 
the safety of flight 

♦> Develop model update techniques based on design optimization to improve analysis/test correlation 



Flutter Analysis Procedure @ NASA Dry den 


□ Flutter Analysis 


❖ Uncertainties in the structural dynamic model are minimized through the use of 
“model tuning technique” 

❖ Based on analytical modes 

Measure Weight, C.G., 
Moment of inertia, 

GVT: <F g ,co g 


□ Validate Structural Dynamic Finite Element Model using Test Data and Update if needed 

❖ Use MDAO (Multidisciplinary Design, Analysis, and Optimization) tool with 
Model Tuning Capability or Standalone Model Tuning Code 

> Model tuning is based on optimization. 

V Design Variables 

• Structural sizing information: Thickness, cross sectional area, 
area moment of inertia, etc. 

• Point properties: lumped mass, spring constant, etc. 

• Material properties: density, Young’s modulus, etc. 

V Constraints 


Structural Dynamic 
Finite Element Model 

I 


MDAO tool with Model 
Tuning Capability 


± 

Validated 

Finite Element Model 

* 


Create Unsteady 
Aerodynamic Model 

* 


Perform Flutter Analysis 
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Test Article: Exploded View of ATW 2 





Structural Dynamic Finite Element Model 


□ Based on MSC/NASTRAN code 

❖ 265 nodes 

❖ Use 10 modes for the flutter analysis 




Structural Dynamics Group 


Chan-gi Pak-5 


Unsteady Aerodynamic Model 


□ Based on ZAERO code 

♦> 416 elements 

Select 16 reduced frequencies between 0 & 1 
Mach = .60, .75, .82, and .95 
Linear Theory 

Use Matched Flutter Analysis 


❖ 


❖ 


❖ 



B 



Mode 1: 17.60 Hz 


V 


\ .'■■'V \ S \. \ \ \. \ 

. \ \ \ \ \ \ \ \ 


Splined Mode Shape 


V \ \ \ \ \ 








\ \ \ \ 




\ \ \ \ \ \ 


\ 


Input Mode Shape 


Structural Dynamic 
Finite Element Model 

Structural Dynamics Group 


V V W WW Aerodynamic 

Model 



Unsteady 


First Bending 



Mode 2: 23.26 Hz 


Nodal Line 



Input Mode Shape 


Structural Dynamic 
Finite Element Model 
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Mode 3: 93.99 Hz 



Structural Dynamic 
Finite Element Model 



Structural Dynamics Group 


Chan-gi Pak-9 


Splined Aerodynamic Loads at Mach = 0.82 


' Ny v \\\\\ \\ \ Unsteady 

\ \ \ \ \ \ \, \ \ \ Aerodynamic 

' ■ X \ ■ - • \ \ \ - \ Model 


Input: Cp Distribution 


Output: Splined Loads 




Structural Dynamic 
Finite Element Model 
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Frequency (Hz) Damping 


V-g and V-oo Curves at Mach = 0.82 Before Model Tuning 



Speed (Keas) 



3% 


Damping 


Mode 

Frequency 

Modal Participation Factor 

1 

17.60 Hz 

75.0 % 

96.6 % 

2 

23.26 Hz 

16.8 % 

3 

93.99 Hz 

4.8 % 

4 

135.4 Hz 

0.0 % 

3.4 % 

5 

163.1 Hz 

2.6 % 

6 

174.5 Hz 

0.0 % 

7 

257.5 Hz 

0.5 % 

8 

391.6 Hz 

0.0 % 

9 

394.3 Hz 

0.1 % 

10 

445.6 Hz 

0.3 % 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

407.4 Keas 

22.86 Hz 

15010 ft 
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Speed (Keas) 
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Summary of the Modal Participation Factors 


Mode 

Frequency 

Modal Participation Factor 

Mach 

= 0.60 

Mach 

= 0.75 

Mach 

= 0.82 

Mach 

= 0.95 

1 

17.60 Hz 

68.1 % 

95.5 % 

72.9 % 

96.2 % 

75.0 % 

96.6 % 

79.7 % 

97.6 % 

2 

23.26 Hz 

22.2 % 

18.3 % 

16.8 % 

13.6 % 

3 

93.99 Hz 

5.2 % 

5.0 % 

4.8 % 

4.3 % 

4 

135.4 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


5 

163.1 Hz 

3.3 % 


2.9 % 


2.6 % 


1.9% 


6 

174.5 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


7 

257.5 Hz 

0.7 % 

4.5 % 

0.6 % 

3.8 % 

0.5 % 

3.4 % 

0.3 % 

2.4 % 

8 

391.6 Hz 

0.0 % 


0.0 % 


0.0 % 


0.0 % 


9 

394.3 Hz 

0.1 % 


0.1 % 


0.1 % 


0.0 % 


10 

445.6 Hz 

0.4 % 


0.3 % 


0.3 % 


0.2 % 



□ Participation of the first three modes is a function of Mach number. 

□ In-plane modes do not participate for the first flutter mechanism at all. 


❖ Modes 4, 6, and 8 

□ Primary Modes: Modes 1, 2, and 3 

❖ Frequency error should be less than 3%. 

□ Secondary Modes: Modes 4 through 10 (higher) 

♦♦♦ Frequency error should be less than 10%. 


Structural Dynamics Group 


Chan-gi Pak-12 



Model Correlation Requirements 


□ Everyone believes the test data except for the experimentalist, and no one believes the finite element 
model except for the analyst. 

♦> Some of the discrepancies come from analytical Finite Element modeling uncertainties, noise in the test results, 
and/or inadequate sensor and actuator locations. Not the same orientation for each sensor. 


□ MIL-STD-1540C Section 6.2.10 

❖ Test Requirements for Launch, Upper-Stage, & Space Vehicles 

❖ Less than 3% and 10% frequency errors for the primary and secondary modes, respectively 

❖ Less than W% off-diagonal terms in orthonormalized mass matrix 


□ AFFTC-TIH-90-00 1 (Structures Flight Test Handbook) 

❖ If measured mode shapes are going to be associated with a finite element model of the structure, it will probably 
need_ to be adjusted, to match the lumped mass modeling of the analysis. 

♦> Based on the measured mode shape matrix <3> and the analytical mass matrix M, the following operation is 
performed: 

<d t m<d 

❖ The results is near diagonalization of the resulting matrix with values close to 1 on the diagonal and values 
close to zero in the off-diagonal terms. Experimental reality dictates that the data will not produce exact unity or 
null values, so 10 percent of these targets are accepted as good orthogonality and the data can be confidently 
correlated with the finite element model. 
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Model Tuning Procedure 


□ Optimization Problem Statement 


□ 


Minimize J 


❖ 


Such that 


J, 


< £, 


A: = 1 ... 1 3 & k*i 


Step 1 : Improve Rigid Body Mass Properties 

❖ Errors in Total Mass 

❖ Errors in CG Locations 


Infeasible 


❖ 


Errors in Mass Moment of Inertias 


Mass Properties 


Total Mass 


CG Locations 


Mass Moment of Inertias 


Objective Functions & Constraints 


J, = (W-W G ) 2 /W G 2 


J 2 = (X-X G ) 2 /X 2 


J, = (Y-Y G ) 2 /Y G 2 


J 4 = (Z-Z g ) 2 /Z 0 2 



Measure Weight, C.G., 
& Moment of inertia 
GVT 

<3>g,cog 


x = (I 


T F/I 2 

A wn ) ' 


XX XXG 


XXG 


= a 


T WT 2 

) / At 


YY A YYG 


L YYG 


x = (I 


i wi 2 

1 77 P_ 7 / Ay 


ZZ ZZG 


ZZG 


A. = (i 


i wi 2 

rvvn y / At 


XY A XYG 


XYG 


*^9 (^YZ - ^YZg) 2 ^YZG 2 


*Ao (^ZX“^ZXg) 2 ^ZXG 2 


SD FE Model 


c 


T 


Start 

"T 


3 


Optimization Step 1 : 
Update Mass Properties 

W, X CG , Y cg , Z cg , Ixx, 
Iyy, Izz, Ixy, Iyz, & Izx 


Optimization Step 2: 
Improve Orthogonality 
of Mass Matrix M 


Optimization Step 3: 
Update Frequencies 
& Mode Shapes 


c 


T 


End 


T 


3 


Aero Model A 

\ 


Flutter Analysis 
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Model Tuning Procedure (continued) 


□ 


Step 2: Improve Mass Matrix 

❖ Off-diagonal terms of Orthonormalized Mass Matrix: M = <E> G T T T M TO G 

n 


J n = 




i= lj=l,i>y 
Guyan reduction 



/ 


M 

M 1 


K 

K 1 

II 

II 

F-h 

— K K 

^ ss ^ sm 

M = 

mm 

M sm 

ms 

Ks 

K = 

mm 
v sm 

ms 

K ss 


Measure Weight, C.G., 
& Moment of inertia 
GVT 

<t> G ,co G 


Improved reduction system 


T = T = 

1 1 IRS 


I 


-K;X m+ (K;lM ss -ft/xwx 

M g =T t g MT g K g =T t g KT g 


System Equivalent Reduction Expansion Process (SEREP) 

r 1 


T = T, 


SEREP 


<d = 


q> kp r q) T <p r 

m V mm/ m 


O lo 7 O ) o 

v \ mm/ m 


o 

n 

<x> 
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FEM 


C 


SD FE Model 


Start 

~r 




Optimization Step 1 : 
Update Mass Properties 

W, X CG , Y cg , Z cg , Ixx, 
Iyy, Izz, Ixy, Iyz, & Izx 


Optimization Step 2: 
Improve Orthogonality 
of Mass Matrix M 


Optimization Step 3: 
Update Frequencies 
& Mode Shapes 



T 


End 


T 


3 
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i 
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Model Tuning Procedure (continued) 
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Test Setup: #1 GVT (Strong Back Mounting) 



Excitation point 




26 27 28 29 3C 


^ \ 

15 16 17 18 19 20 21 22 23 24 25 I 

Excitation point 


PONTOS photogrammetry 
optical measuring system 
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Test Setup: #2 GVT (FTF Mounting) 



Impulse Hammer 


1 
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Flight Test Fixture 


ATW2 





Test Setup: #3 GVT (FTF Mounting under F-15) 



Structural Dynamics Group 


Chan-gi Pak-19 



Sensors 



Teardrop accelerometers 
(2: Epoxied onto Surface) 
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Results (Frequency Comparisons) 


Mode 

GVT (Hz) 

Before Tuning 

After Tuning (Target #2 GVT data) 

MIL- 

STD 

(%) 

#1: 

Strong 

Back 

#2: 

FTF 

#3: 

FTF & 
F15B 

Freq 

(Hz) 

Error (%) 

Freq 

(Hz) 

Error (%) 

Wrt #1 

Wrt #2 

Wrt #3 

Wrt #1 

Wrt #2 

Wrt #3 

1 

17.24 

17.45 

17.42 

17.60 

2.09 

0.86 

1.03 

17.45 

1.22 

0.00 

0.17 

3 

2 

44.10 

43.72 

43.73 

23.26 

- 47.3 

- 46.8 

- 46.8 

43.48 

- 1.41 

- 0.55 

- 0.57 

3 

3 

84.00 

83.66 

84.14 

93.99 

11.9 

12.4 

11.7 

82.98 

-1.21 

-0.81 

-1.38 

3 

4 

N/A 

N/A 

N/A 

135.4 

N/A 

N/A 

N/A 

133.6 

N/A 

N/A 

N/A 


5 

N/A 

142.3 

143.0 

163.1 

N/A 

14.6 

14.1 

153.8 

N/A 

8.08 

7.55 

10 
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Results (Total Weight, Orthogonality, & MAC) 
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Frequency (Hz) Damping 


V-g and V-co Curves at Mach = 0.82 After Model Tuning 



Mode 

Frequency 

Modal Participation Factor 

1 

17.45 Hz 

10.7 % 


2 

43.72 Hz 

87.2% 

99.5 % 

3 

83.66 Hz 

1.6 % 



Speed (Keas) 
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200 300 400 500 600 700 


Flutter Mode 

Speed 

Frequency 

Altitude 

1 

341.5 Keas 

34.59 Hz 

23475 ft 
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Flutter Results Before & After Model Tuning 


Mode 

Frequency 

Modal Participation Factors before Model Tuning 

Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

1 

17.60 Hz 

68.1 % 

95.5 % 

72.9 % 

96.2 % 

75.0 % 

96.6 % 

79.7 % 

97.6 % 

2 

23.26 Hz 

22.2 % 

18.3 % 

16.8 % 

13.6 % 

3 

93.99 Hz 

5.2 % 

5.0 % 

4.8 % 

4.3 % 

Mode 

Frequency 

Modal Participation Factors after Model Tuning 

Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

1 

17.45 Hz 

5.0 % 

99.7 % 

8.2 % 

99.6 % 

10.7 % 

99.5 % 

22.6 % 

96.4 % 

2 

43.72 Hz 

93.8 % 

90.0 % 

87.2 % 

71.4 % 

3 

83.66 Hz 

0.9 % 

1.4% 

1.6% 

2.4 % 




Mach = 0.60 

Mach = 0.75 

Mach = 0.82 

Mach = 0.95 

Before 

Tuning 

Speed 

453.0 Keas 

421.5 Keas 

407.4 Keas 

377.9 Keas 

Frequency 

23.18 Hz 

22.97 Hz 

22.86 Hz 

22.53 Hz 

Altitude 

-7501 ft 

8751 ft 

15010 ft 

25590 ft 

After 

Tuning 

Speed 

337.9 Keas 

340.5 Keas 

341.5 Keas 

344.7 Keas 

Frequency 

35.96 Hz 

35.11 Hz 

34.59 Hz 

32.91 Hz 

Altitude 

8642 ft 

19400 ft 

23475 ft 

29700 ft 

Speed Error 

33.8 % 

23.8 % 

19.3% 

9.6 % 
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Altitude (ft) 


Flutter Boundary vs. Flight Envelope 


50000 


40000 


30000 


20000 


10000 


0 


-10000 
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: Flutter Boundary Before Model Tuning (3% structural damping) 

: Flutter Boundary After Model Tuning (3% structural damping; reference 2) 

: Flutter Boundary After Model Tuning (3% structural damping; FEM based on #2 GVT data) 
: Flutter Boundary After Model Tuning (measured damping; FEM based on #2 GVT data) 
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: Mother Ship Envelope 


: ATW2 Measured Flutter Boundary 
ATW2 Measured Flutter Boundary/ 1.15 


0.5 


Mach 



1.5 
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Conclusion 


□ Model tuning based on GVT data is needed to minimize uncertainties in the structural 
dynamic model and to increase the safety of flight. 


□ After model tuning (for ATW 2 case) 

♦♦♦ Maximum of 47% primary frequency error in second mode becomes 0.55%. 

> Should be less than 3% 

♦♦♦ Maximum of 66% error in off-diagonal terms of orthonormalized mass matrix becomes 6.2 %. 

> Should be less than 1 0% 


♦♦♦ Maximum error in flutter speed was 33.8% (Mach = 0.60) 


> 

> 


Modal participation of the second mode was more dominant than the first mode after model tuning. 


Modal participation of the second mode vs. Flutter speed error 



Mach = 0.60: 

93.8% 

33.8% 


Mach = 0.75: 

90.0% 

23.8% 


Mach = 0.82: 

87.2% 

19.3% 


Mach = 0.95: 

71.4% 

9.6% 
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Conclusion (continued) 


Flutter Boundaries 

Flutter Speed Differences 

Measured/ 1 . 1 5 = V d 

0% 

Measured = 1 . 1 5 V d 

15% 

FEM; after model tuning; based on #2 GVT data; with measured damping 

37% 

FEM; after model tuning; based on #2 GVT data; with 3% structural damping 

41% 

FEM; after model tuning; based on #1 GVT data; with 3% structural damping 

52% 

FEM; before model tuning; with 3% structural damping 

74% 


Validated Structural 
Dynamic Model 

Validated Unsteady 
Aerodynamic Model 

Recommended* 
Flutter Margins 

ATW2 Case 

Yes 

Yes 

15% 


Yes 

No 

49% 

40 - 50 % 

No 

No 

54% 

74% 


*: JSSG-2006 Guidelines for Flutter Speed Clearance; Faustino Zapata, AFDC May 22-23, 2008 
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